Scaling properties of scale-free evolving networks: Continuous approach 
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Scaling behavior of scale-free evolving networks arising in communications, citations, collabo- 
rations, etc. areas is studied. We derive universal scaling relations describing properties of such 
networks and indicate limits of their validity. We show that main properties of scale-free evolving 
networks may be described in frames of a simple continuous approach. The simplest models of 
networks, which growth is determined by a mechanism of preferential linking, are used. We con- 
sider different forms of this preference and demonstrate that the range of types of preference linking 
producing scale-free networks is wide. We obtain also scaling relations for networks with nonlinear, 
accelerating growth and describe temporal evolution of arising distributions. Size-effects - cut-offs 
of these distributions - implement restrictions for observation of power-law dependences. The main 
characteristic of interest is so-called degree distribution, i.e., distribution of a number of connec- 
tions of nodes. A scaling form of the distribution of links between pairs of individual nodes for the 
growing network of citations is also studied. We describe effects that produce difference of nodes. 
"Aging" of nodes changes exponents of distributions. Appearence of a single "strong" node changes 
dramatically the degree distribution of a network. If its strength exceeds some threshold value, the 
strong node captures a finite part of all links of a network. We show that permanent random damage 
of a growing scale-free network - permanent deleting of some links - change radically values of the 
scaling exponents. We describe the arising rich phase diagram. Results of other types of permanent 
damage are described. 
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I. INTRODUCTION 

One of the most impressive recent discoveries in the 
field of the network evolution is an observation that a 
number of large growing networks are scale-free, that is 
distributions of numbers of connections of their nodes 
are of a power-law form [l]-[|] . To the class of scale- free 
networks belongs huge communications networks (World 
Wide Web and Internet networks of citations in 

literature H|l(| , collaboration networks Jll],[l^] , some bi- 
ological networks (nets of metabolic reactions |L3]]), etc. 
In fact, an incredible progress of science and information 
technology itself has produced networks, large enough 
(World Wide Web has about 10 9 nodes) to get reliable 
data. 

Nevertheless, the experimental data are certainly not 
excellent. Indeed, the task of observation of such power 
laws is not easy. Arising distributions are very sensible 
to size-effects. Unfortunately, few really large growing 
networks are known yet, many of obtained data for not 
so large nets are very preliminary. It is a great tempta- 
tion to describe any decreasing experimental curve in a 
log-log plot by a linear fit. Nearly all the observations 
are for the only quantity - degree distribution of nodes. 
(Following definitions of mathematicians and computer 
scientists, we call here a total number of connections of a 
node its degree |l4| ] . If the links are connected, a number 
of incoming links of a node is called in- degree of it, and 
a number of outgoing links - out- degree.) Information 
about others characteristics is much poorer. 



Therefore, only a few scale-free growing networks have 
been observed. Why are these observations so impor- 
tant? Why the interest in these networks is so great? (In 
one of the last issues of Physical Review Letters, three 
papers, one after the other, were devoted to scale-free 
networks [[l5-17|.) Of course, the reason is not only 
power-law dependences of distributions itself but a va- 
riety of intriguing properties of scale- free networks which 
explain their existence in Nature. Degree distribution is 
a simple but very important characteristic of a network. 
In particular, it was shown that networks with power- 
law degree distributions are extremely resilient against 
random breakdowns if an exponent 7 of their degree dis- 
tribution does not exceed 3. This property is vitally 
important for communications and biological networks. 
One has also point out the following circumstance. A 
power-law dependence of degree distributions indicates 
scale-invariance of characteristics of networks and their 
hierarchically organized structure. Just investigation of 
scaling properties of scale-free networks and connections 
between them is a topic of the present paper. 

We study the formulated problem applying two dif- 
ferent approaches in arbitrary order. First, we obtain 
universal scaling relations using general considerations. 
Second, we demonstrate features of scaling behavior of 
networks using minimal models providing the effect be- 
ing studied. At the moment, the only known mechanism 
providing scale-free networks is preferential linking new 
links are preferentially attached to nodes with large num- 
ber of connections (large degree) ]5|, ^8| , ^9| . Probability 
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that a new link turns to be attached to a node is a func- 
tion of its degree. Related ideas was discussed long ago 
by Simon pfl]. All these models belong to a class of 
stochastic multiplicative processes [Ell. 

All the used models can be described by linear mas- 
ter equations Jl7j which are solvable exactly in several 
particular cases. In this pa per, we treat them using a 
trivial continuous approach |!8|Jl9|Jl7[| that, as we show, 
describes main features of the networks and produces ex- 
act values of the exponents. Demonstration of possibili- 
ties of the continuous approach is one of the aims of our 
paper. 

In our paper, we use several different models of evolv- 
ing networks (with undirected and directed links), but a 
skeleton of all them is the same (see Fig. |l|). At each 
increment of time, a new node is added. Therefore, the 
number of the last node, t, may be called "time" . Each 
node is labeled by the time of its birth, s = 0, 1, 2, . . . < t. 
Together with a new node, several new links is appeared 
in the network. A part of these links is distributed be- 
tween nodes preferentially according to their degrees (we 
notate it k) or in-degrees (= q). (Out-degree distribu- 
tions is discussed in only one place). In general consid- 
erations, k and q will be equivalent if no special notion 
will be made. 

The following "one-node" characteristics may be in- 
troduced. p(k, s, t) is the probability that the node s has 
degree k at time t. P(k, t) is total degree distribution 
of a network in time t. P(k) — P(k, t — » oo) is station- 
ary distribution, 7 is the corresponding scaling exponent, 
P{k) oc A; -7 for large k. k(s, t) is the average degree of the 
node s at time t. (3 is its scaling exponent, k(s, t) oc s^ 13 . 
The corresponding definitions for in-degree distributions 
are similar. 

In Sec. |TJ we introduce the continuous approach. Con- 
sidering a simple example we explain reasons to use this 
approach and derive useful relations. 

In Sec. Ill, for linearly growing scale- free networks (to- 
tal number of links is proportional to a total number of 
nodes with time-independent coefficient), using general 
considerations, we obtain scaling forms of involved quan- 
tities and universal relations between scaling exponents. 

Sec. |y] is devoted to study of different types of pref- 
erential linking. Using the continuous approach, we de- 
scribe (i) the simplest network in which new links are 
attached without any preference (produces exponential 
distributions), (ii) linear type of preferential linking (pro- 
duces scale- free networks with 2 < 7 < 00), (iii) a mix- 
ture of preferential and random linking (also produces 
scale- free networks with 2 < 7 < 00), and (iv) more gen- 
eral cases also producing scale-free networks. In fact, in 
Sees. [jv|, [v], and VI different realizations of a linear pref- 
erence function, G(s, t)k+A(s, t), according to which new 
links are distributed among nodes, are studied. In Sec. 
[V|, "strength" of a node, G(s, t), depends only on its age, 
G(s,t) oc (t — s)~ s . This changes scaling exponents. In 
Sec. [Vl| strength of nodes depends only on dates of their 
birth, s. In this case, the network may exibit an intrigu- 



ing phenomenon of capture of a finite fraction of all links 
by the strongest node - giant capture, - that was quite re- 
cently reported p2fl . Nevertheless, we demonstrate that 
the network stays to be scale-free. 

Note that, in most of considered here models, new links 



may also connect old nodes. Nevertheless, in Sec. VII 



we study a specific growing network in which new con- 
nections are possible only between a new node and old 
ones, that, in particular, may be used for description of 
networks of scientific citations. In this simple situation, 
in the continuous approach, it is possible to describe dis- 
tributions of links between pairs of individual nodes. 
Nonlinearly growing networks are considered in Sec. 



VIII. For this, more general case, we obtain scaling rela- 
tions, and describe non-stationary distributions, P(q,t), 
and their cut-offs impeding observation of power-law de- 
pendences for finite-size networks. 

In Sec. IX, we study influence of permanent damage on 
the scaling characteristics of growing networks. We show 
that such type of the damage produces much stronger 
effect on the network than a previously studied instant 
damage ^g|Mp||. 

High quality of the used continuous approach is dis- 
cussed in Sec. [X| using already known exact results. 



II. REASONS FOR THE CONTINUOUS 
APPROACH 

We start from one of the most simple models of growing 
networks with preferential linking proposed by Barabasi 
and Albert H which belongs to classes of more general 
models wich were solved exactly afterwards |l^Jl^]. At 
each increment of time, a node is appeared. It connects 
to one of old nodes chosen with probability proportional 
to degree of this old node, i.e., to a total number of its 
connections. Using introducing in Sec. [| notations, it is 
possible to write immediately the following master equa- 
tion for degree probabilities of individual nodes (see jlTj ) , 

p{k, s,t+l) = ^-P(k - 1, s, t)+(l- p(k, s, t) 

(1) 

with, e.g., t = 1,2,3,... and s — 0, 1,2, i.e., at 

time t = 1, a pair of nodes, s — 0, 1, connected by 
a link is present. Therefore, the initial condition is 
p{k,s = 0, l,t = 1) = Sk,i and p(k,t,t) — Sk,i is the 
boundary condition. Eq. ([!]) may be rewritten in the 
form, 

2t\p(k,s,t+l) -p(k,s,t)) = 

(k-l)p(k-l,s,t)-kp(k,s,t). (2) 
Passing to the continuous limit in t and k we get 



2t dp(k,s,t) d[kp(k,s,t)] = Q 



dt 



dk 



(3) 
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and 



d[kp(k, s, t)\ d[kp(k,s,t)} = Q 



dlnVi 



dink 



(4) 



The solution of Eq. (Q) is kp(k,s,t) = 5(\nk - 
In \ftjs~ + const) . The boundary condition is fulfilled if 
it is of the following form, 



p(k, s, t) = S(k - 



(5) 



Therefore, we see that the transition to the continuous 
limit in the master equation leads to the ^-function form 
of degree distributions of individual nodes. 

The main quantity of interest is total degree distribu- 
tion of all the network, 



The meaning of Eq. (|T(]) is quite obvious - each 
new link is distributeded homogeneously among all nodes 
with account of the particular form of preference. At 
first, this approach has been called "mean field" PJl^] 
but afterwards it was shown that it is simply continuous 
approximation |l^ , |l7|] . 

One should note that we pass to the continuous limit 
for both variables - k and t. The continuous limit for t is 
not dangerous if we consider large enough networks and 
do not study some peculiarities of distributions related 
with particular form of initial conditions (for a very sim- 
ple model of a scale free network, an exact solution was 
found for all values of k and t |26|). The continuous in k 
approximation needs more care, so we discuss its quality 
along all the paper. 
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III. SCALING RELATIONS FOR SCALE-FREE 
NETWORKS 



In the continuous approximation, stationary degree dis- 
tribution is of the form, 



P(k) = P{k,t^> oo) = lim - [ dsp(k, 



8,t). 



(7) 



Inserting the obtained expression for p(k,s,t), Eq. (|^), 
into Eq. (|^), one gets the continuous approximation re- 
sult for this model, P(k) = 2/fc 3 |@. 

Another way to obtain this expression is to derive 
an equation for the total probability, P(fc,f). Applying 
X)!=o t° both sides of Eq. (]]]) and passing to the t — > oo 
limit, one gets 

P{k) + \[kP{k) - (k- l)P(k - 1)] = S kA . (8) 

In the continuous k limit, this equation is of the form 
P(k) + (l/2)d[kP(k)]/dk = 0. Its solution is P(k) oc fc~ 3 . 

Often, it is more simple to proceed in a bit different 
way. Let us introduce average connectivity of an individ- 
ual node: 



k(s,t) = y2kp(k,s,t) = I 
k=i Ja 



dk kp(k, s, t) . 



(9) 



Applying f °° dkk to Eq. (Q) and integrating its right 
part by parts, we get 



dk(s,t) 
dt 



k(s,t) 



J du k(u, t) 



(10) 



Here, we used that, in the particular case of this 
model, the total number of links in the network, 
Jq* ds dk kp(k, s, t), equals It. Boundary condition for 
Eq. (|l(| is k(t, t) = 1. It is usefull to use equations simi- 
lar to Eq. ([h]) but not to study each time corresponding 
master equations. 



In the continuous approach, knowledge of average de- 
gree of nodes, k(s, t), lets to obtain the total degree distri- 
bution, P(k, t). (Now we do not restrict ourselves by any 
particular model and our general results are also valid for 
in-degree, q, (out-degree) distributions. Nevertheless, in 
any case, q(s,t) is the solution of an equation similar to 
Eq. ©.) 

Indeed, one can write 
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P(k,t) = - / dsS(k- k(s,t)) 
t Jo 



1 fdk(s,t) 



t V ds 



[s = s(k,t)}, 



(11) 



where s(fc,t) is a solution of the equation, k = fc(s,t). 

Now, we can easily connect involved quantities. As- 
suming that P{k) and k(s, t) exibit scaling behavior, that 
is P(k) oc fc~ 7 for large k and k(s, t) oc s _/3 for 1 <C s <C t, 
one gets s oc k" 1 ^ and fc~ 7 oc ds/dk oc /c~ 1 ~ 1 /' 3 . There- 
fore, 7 = 1 + 1//3, and we obtain the following scaling 
relation: 



0(7-l) = l 



(12) 



Let us show that the relation, Eq. (|l2|), is univer- 
sal. Here, we can proceed using general considerations. 
In the present section, we study only linearly growing 
networks (input flow of new links does not depend on 
time) which produce stationary degree distributions at 
long times. More complex cases will be considered in 
Sec. |VIIIA . 



Hence, let P(k) be stationary, then it follows from 
Eq. (Q) that p(k, s, t) has to be of the form, p(k, s, t) = 
p(k,s/t). The normalization condition is 



dkp(k, s, t) = 1 



(13) 
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so J °° dk p(k, x) = 1. From this equation, it follows that 
p{k,x) = g(x)f(kg(x)), where g{x) and f(x) are arbi- 
trary functions. From the definition of k(s,t), Eq. (|^), 
using its scaling behavior, we get J Q dkkp(k,x) oc x~P. 
Substituting p{k 1 x) into this relation, one obtains that 
g(x) oc x 13 . Of course, without loss of generality, one may 
set g{x) = x@ , so we obtain the following scaling form of 
degree distributions of individual nodes, 



p(k,s,t) = (s/tff(k(s/tf) 



(14) 



Finally, assuming scaling behavior of P(k), i.e., 
J °° dx p(k, x) oc fc -7 , and using Eq. (|T4|), we obtain 
7 = 1 + 1//3, i.e., the relation, Eq. ( [l2| ) is universal. Here 
we used quick convergence of p(k, x) at large x (it follows 
from our exact results JTtJ). One should note that, during 
this derivation, we did not use any approximations. 



IV. TYPES OF PREFERENTIAL LINKING 

A. Absence of preference 

We need an example of a non-scale-free network for 
comparison, so we start with the simplest growing net- 
work with random attaching of new links. Let again, as in 
Sec. [□], at each increment of time, a new node be added 
to the network. Now it connects with a randomly chosen 
(i.e., without any preference) old node. The growth be- 
gins from the configuration consisting of two connected 
nodes at time t = 1. This model is even simpler than 
the well-known Erdos-Renyi model [^,^8) , and its exact 
solution is trivial. The master equation describing evo- 
lution of degree distributions of individual nodes looks 
as 

p(k,s,t+ 1) = —^—p(k - l,s,t) + ( 1- — !— ) p(k,s,t) . 



t+V 



t + l 



(15) 



p(k,s = 0, l,t = 1) = S(k,t,t> 1) = 4,i- Applying 



J2 s =o t° both sides of Eq. ( [i5| ) and using the definition 
of total degree distribution, Eq. (||), we get 

(t + l)P(k, t + l)- tP{k, t) = P(k - 1, t) - P(k, t) + 4,i 

(16) 

The corresponding stationary equation, 

2P{k) - P(k - 1) = 4,1 , (17) 

has a solution of an exponential form, P(k) = 2~ k . 
Therefore, networks of such type often are called "ex- 
ponential" . 

From Eq. ([l5|), one may also find the degree distribu- 
tion of individual nodes, p(k, s,t), for large s and t and 
fixed s/t, 



p(k,s,t) 



1 



t{k + l)\ 



In* 



fe+i 



(18) 



One sees that this function decreases rapidly at large val- 
ues of degree, k. 

Thus, degree distributions of networks growing with- 
out using a mechanism of preferential linking differs strik- 
ingly from distributions of scale-free networks described 
in previous section. Note that the singularity of p(k, s, t), 
Eq. (|l8|), at s/t — > is much weaker than for scale- free 
networks. 



B. Linear preference 

Let us demonstrate the continuous approach in details 
using a more general model than in Sec. O which pro- 
duce a wide range of 7 exponent values JlTj (see Fig. |l , a) . 
Let us consider a network with directed links. Here, we 
study distributions of numbers of incoming links of nodes 
{in-degree, q). 

At each time step a new node is added. It has n in- 
coming links. These links go out from arbitrary nodes 
or even from some external source. Simultaneously, m 
extra links are distributed with preference. This means 
that again they go out from non-specified nodes or from 
an external source but a target end of each of them is 
attached to a node chosen preferentially: probability to 
chose some particular node is proportional to a function 
of its in-degree, we call it a preference function. In the 
simple model that we consider in the present section, this 
probability is proportional to q+ A, where A is a constant 
which we call additional attractiveness. We shall see that 
its reasonable values are A > — n. 

In the continuous approach, we may assume that m 
and n are not necessary integer numbers but any posi- 
tive ones. We do not worry about source ends of links 
because study here only in-degree distributions. 

Then the equation for average in-degree is of the form, 



dg(s,t) 
dt 



q(s,t)+A 
f* du[q(u,t) + A] 



(19) 



with the initial condition, q(0, 0) — 0, and the boundary 
one, q(t, t) = n. 



Applying J Q ds to Eq. (19) we obtain 

d [* f t d_ 
— J dsq(s,t) = J ds —q(s,t) + q(t,t) , 



so 



ds q(s, t) = (to + n)t . 



(20) 



(21) 



This relation is quite natural - the total in-degree of the 
network is equal to the total number of created links, and 
we see that Eq. ( |Ts| ) is consistent. Therefore, Eq. ( |l9| ) 
takes the form, 
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dq(s,t) = _ 
dt to 

Its general solution is 



m q(s, t) + A 
-n + A t 



l{s, t)+A = c{s)t m/{m+n+A) . 



(22) 



(23) 



where C(s) is arbitrary function of s. Accounting for the 
boundary condition, q(t, t) = n, one has 



q(s,t) + A=(n + A)^ 

Hence, the scaling exponents are 

m 



-raj (m+n+A) 



P = 



m + n + A 



and 



7 = 2- 



A 



(24) 



(25) 



(26) 



One sees that, for n + A > 0, the exponent 7 is in 
the range 2 < 7 < 00 while f3 belongs to the interval, 
< (3 < 1. We see that n plays the same role as A. 

If we set n = and demand that all new links have 
to go out from new nodes, we get a network of citations 
(see Fig. [l], b). Note that, when A = m and n = 0, we 
get the particular case of the Barabasi- Albert's model |jj 
in which each new link is connected with a new node. 
Indeed, in this case, degree and in-degree of nodes are 
coupled rigidly, k — q + to, and the obtained results are 
valid for the degree distributions. If, in addition, we set 
m = 1, we obtain the model considered in Sec. O. 



C. Mixture of preferential and random linking 

Now we can consider a bit more complex model. We 
will demonstrate that scale-free networks may be ob- 
tained even without "pure" preference linking. We dis- 
cuss the model introduced in the previous section but 
with one new element. Let, in addition, at each time 
step, n r links be distributed between nodes randomly, 
without any preference. Again, we study in-degree dis- 
tributions, so these links may go out from anywhere but 
their target ends are attached to randomly chosen nodes. 
Then, 



dqjs, t) 
dt 



n. 

— + TO 



q(s,t)+A 



t f*du[q(u,t) + A] 



(27) 



Again, q(0, 0) = and q(t, t) = n. In this case, 
/J dsq(s, t) = (n r + m + n)t, so Eq. ( ^7|) takes the form: 

dq(s, t) to 



dt 



n r + m + n + A 



q(s,t) + A+ —(n r + to + n + A) (28) 
to J 



Its solution, with account for the boundary condition is 



q(s, t) + A + — (n r + m + n + A) = 

TO 

n r +m, / s ^-m/(n r +m+n+A) 

— (n + nr + A) - . 29 

TO V t J 



Here, to, n, and n r are positive numbers, and n r +n+A > 
0. 

Therefore, we get the scale-free network with the ex- 
ponents, (3 = to/(to + n r + n + A) and 



7 = 2- 



n r + n + A 



(30) 



Thus, additional fraction of randomly distributed links 
does not delete a power-law dependence of distributions 
but only increase 7. 

In two computer-science papers [^9],U a similar model 
of a network with directed links was considered. At each 
time step, a new node is added to the net. It has one out- 
going link. The other end of this link connects with one 
of old nodes by the following rule, (i) With probability 
p it connects with a random node, (ii) With probability 
1 — p it finds a random node and connects with its sole 
target neighbor node. (In this particular case, it is the 
same as to chose random link and to connect with its 
target end.) 

One can see that this model corresponds to the par- 
ticular case, n = 0,n r = p,m = 1 — p, A — 0, of the 
network that we consider in the present Section. From 
Eq. ©, we get 7 = 2 + p/(l - p) = 1 + 1/(1 - p). 
(Note that, in these papers Refs. [[29P1, unproper result 
figurates. The necessary additional unit is absent there. 
Indeed, if one of the factors (random linking) produces 
7 = 00, another one (preferential linking) - 2 < 7 < 00, 
then their interplay can not produce 7 < 2.) 

It is worth while to emphasize here a close connection 
of the models that we consider with the well-known Si- 
mon's model pc| discussed in relation with networks in 
p0{ . In the Simon's model, if one uses the language of 
growing networks with directed links, at each time step, 
a new link appears. Since we discuss here only in-degree, 
it is again not important, from were it goes out. With 
some fixed probability, a new node is created and the 
target end of the link is attached to it. With an adjusted 
probability, the target end of the link is attached to the 
target end of a randomly chosen old nodes. Then only 
one point in the Simon's model, differs from the models 
that we consider in the present paper: at each increment 
of time, a link is added but not the node as in our case. 
Of course, this does not influence results for large net- 
works and scaling exponents. 
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D. WWW exponents 



The previous result for the 7 exponent, Eq. (p0|), lets 
us to obtain the most rude estimates for the exponents 
of the World Wide Web |}0 31 1 which are measured with 
sufficient precis ion f or comparison. Let us apply the 
model of Sec. 



[V C to the growth of the World Wide 



Web. This means that we have to assume that each time 
when a new node of the Web appears, in average, the 
same number of new links arises between its nodes. We 
neglect many very important factors including eliminat- 
ing of some nodes and links during the growth, etc. 

We do not know values of any of the quantities in the 
left part of Eq. ([50|). The constant A may takes any 
values between — (n r + n) and infinity, the number of 
the randomly distributed links, n r , in principle, may be 
not small (there exist many persons making their refer- 
ences practically at random), n is not fixed. From the 
experimental data [Q, we know more or less the sum, 
m + n + n r ~ 10 3> 1 (between 7 and 10, more precisely), 
and that is all. 

The only thing we can do, it is to fix the scales of the 
quantities. The natural characteristic values for n r +n+A 
in Eq. @ are (a) 0, (b) 1, (c) m > 1, and (d) infinity. 
In the first case, each node have zero initial attractive- 
ness, and all new links are directed to the oldest node, 
7^2. In the last case, there is no preferential linking, 
and the network is not scale- free, 7 — > 00. Let us consider 
the really important cases (b) and (c). 

(b) How do pages appear in the Web? Suppose, you 
want create your personal home page. Of course, first 
you prepare it, put references, etc. But that is only the 
first step. You have to make it accessible in the Web, to 
launch it. You come to your system administrator, he 
put a reference to it (usually one reference) in the home 
page of your institution, and that is more or less all - 
your page is in the World Wide Web. There is another 
way of appearing of new documents in the Web. Imag- 
ine that you already have your personal home page and 
want to launch a new document. The process is even 
simpler than the one described above. You simply insert 
at least one reference to the document into your page, 
and that is enough for the document to be included in 
the World Wide Web. If the process of appearing of each 
document in the Web is as simple as the creating of your 
page - only one reference to the new document (n = 1) 
- and if one forgets about the terms n r and A in Eq. 
(|30|), than, for the 7j„ exponent of the distribution of 
the incoming links (in-degree distribution), we get im- 
mediately the estimation, 7j„ — 2 ~ 1/m ~ 10 _1 . This 
estimation is indeed coincides with experimental value 
ji n - 2 = 0,1 0. (Here, we have introduced the new 
notation, 7j„, because in the present section we discuss 
different distributions.) Therefore, the estimation seems 
to be good. Nevertheless, we should repeat again, that 
this estimation follows only from the fixation of scales 
of the involved quantities. We emphasize that there are 



no any general reasons to set, e.g., A — 0. A lot of real 
processes are not included in this estimation. Aging of 
nodes changes 7 see Sec. [vjand account for dying of 
nodes 7, see |32) and (the half- life of a page in the Web 
is of the order of half a year) changes 7. The ratio be- 
tween the total number of links and the number of nodes 
in the Web is not constant [Q , it increases with time, the 
growth of the Web is nonlinear. This factor also change 
the value of 7, in futu re it m ay become even lower than 
2, see 



and Section VIII 



(c) Above we discussed the distribution of incoming 
links. Eq. (|3(]) may be also applied for the distribution 
of links which come out from documents of the Web, 
since the model of the previous section may be reformu- 
lated for outgoing links of nodes. In this case all the 
quantities in Eq. ( |30| ) takes other values which are again 
unknown. Nevertheless, one may think that the num- 
ber of the links distributing without any preference, n r , 
is not small now. Indeed, even beginners proceed with 
linking of their pages. Hence, n r ~ m - we have no an- 
other available scale, - and -f ou t — 2 ~ m/m ~ 1. Wc 
again can compare this estimation with the experimental 
value, 7out - 2 = 0.7 j§. 



E. Generalized form of preference 

What are the other forms of preference that produce 
scale-free networks? Let us list the main possibilities. In 
the present section and in Sees. [v| and [v| we consider 
only a linear form of a preference function. Nonlinear 
preference functions are discussed in Sec. |x|. 

The reasonable (linear) form of the probability for a 
new link be attached to a node s at time t is p Sl < = 
G Sl tk Sl t+A Si f. The coefficient G St t may be called strength 
of a node (in Refs. ^,||) it is called "fitness") A Sit is ad- 
ditional attractiveness. 

One can consider the following particular cases. 

(i) G = const, A — A(s). In this case, the additional 
attractiveness, A(s) may be treated as ascribed to indi- 
vidual nodes. A possible generalization is to make it a 
random quantity. One can check that the answers do not 
change crucially - one only has to substitute the aver- 
age value, A, instead of A, into previous expressions for 
scaling exponents. 

Note that n and m may be maked also random and 
substituted by n and m in the expressions for the expo- 
nents. 

There is a more interesting possibility - to construct 
a direct generalization of the network considered in Sec. 
IV C. For this, we may ascribe the additional attractive- 



ness not to nodes but to new links and again make it 
random quantity. Therefore, new links play the role of 
fans with different passion for popularity of their idols, 
nodes. This is case (ii), G = const, A — A(t), where A(t) 
is random. 



G 



Let a distribution function of A be P(A). Then, our 
main equation looks as 



dt 



J du[q(u, t) + A] 



(31) 



Initial and boundary conditions are g(0, 0) = and 
q(t,t) =n. Hence, J Q dsq(s,t) — (n + m)t. Therefore, 



dq(s,t) m 
dt = T 



q(s,t) 



dAP(A) 
n + Wi + A 



dAAP(A) 
n + m + A 



Hence we again obtain a scale-free network, 

_ f dAP(A) 
p = m / — — — - 

J n + m + A 



and 
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1 



dAP(A) 
1 + (n + A)/m 



(32) 



(33) 



(34) 



These expressions generalize the corresponding results of 
Sees. [VB| and IV C. One can see that the values of the 
7 exponent are again between 2 and oo, /3 is between 
and 1. 

(iii) A = const. 

The case, G(s,t) oc (s — t)~ A is considered in Sec. [v|, 
it produces 2 < 7 < 00. 

Different versions of the s-dependent strength, 
G(s,t) = G(s), are considered in Sec. 
of the power-law dependence, G(s) 



VI 

-A 



The case 



oc s~" is considered 
in Sec. VIA - three different types of behavior are possi- 
ble: A' < - an exponential network; A' = - scale-free; 
A' > - the oldest node gets a finite part of total degree 
of a network. 

One can consider a network in which several nodes are 
"stronger" then other s. We investigate effects arising in 
this situation in Sec. VI B| . A homogeneous mixture of 



nodes with two different strengths is studied in Sec. VI C 



The case of fluctuating G(s) is considered by Bianconi 
and Barabasi |34j - when the distribution of G is homoge- 
neous, i.e., G is homogeneously distributed between two 
fixed values, the distribution has a logarithmic factor, 
P(q) oc q 1 /hiq, 2 < 7 < 00. For some special forms 
of the distribution of G, strong cooperative effects were 
found recently |22j |. 

(iv) A = const, G — G(t). This case is reduced to the 
case (ii). 



V. EFFECT OF AGING OF NODES 

How does the structure of the network change if one 
introduces aging of sites |19| , i.e., if the probability of 
connection of the new site with some old one is propor- 
tional not only to the connectivity of the old site but 



also to the power of its age, r _A , for example? Here, we 
introduce the aging exponent, A. Such aging is natural 
for networks of scientific citations in which old papers 
usually have low attractiveness. 

The simplest analytical expressions may be obtained 
for the model of Sec. ||: links are undirected, and every 
new node is connected by a single link with some old node 
which is chosen according to the same rule of preference 
as in Sec. y. Fig. || demonstrates that the structure of 
the network depends crucially on the value of the aging 
exponent. When A is large enough (larger than I), the 
network turns to be a linear structure. For low negative 
values of A, i.e., for large —A (attractiveness of docu- 
ments increase with growth of their age) many links are 
attached to the oldest nodes. 

The resulting equation for such a network is of the 
form, 



dk(a,t) k(s,t)(t-s)- 



dt 1 - 

Jduk(u,t)(t - u)- A 




k{t,t) = \. (35) 



Then, f*dsk(s,t) = It. 

We search for the solution of Eq. (^ 
form, 

k(s,t) = n(s/t) , s/t = £, . 



Then Eq. @ becomes 
A dhiK(£) 



in the scaling 



(36) 



k(1) = 1, 



1 I- 1 

-A 



dC«(C)(i-0" 



E/?, 

(37) 



where (3 is a constant which is unknown yet. One can 
understand that it is just the exponent (3 of average con- 
nectivity of individual nodes since, in the scaling region, 
£ 1, Eq. (|37|) provides the power law dependence 

K (0 £ We also get the relation f Q c?C K (C) — 2- Our 
aim is to find (3. For this, we have to find a solution of 
Eq. (|3^ ) containing this unknown parameter (3 and sub- 
stitute it into its definition, the left part of Eq. (|37|), or 
into the relation for total number of links. 
The solution of Eq. @ is 



k(£) = -Bexp 



-0 



(38) 



where B is a constant. The indefinite integral in Eq. (|3q ) 
may be taken: 



I t(l-t) A 

OO 1 

1 ^ + E-tt-tTp a ( a + 1 )---( a + ^ +1 



3=0 



ln£ + A 3 F 2 (l,I,I + A;2,2;£), 



(39) 
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where 3-F2( ,,;,;) is the hypergeometric function. 
Recalling the boundary condition = 1, we find the 
constant B. Thus the solution is 

«(£) = 

_-^[C+V-(l-A)]t-/3 



get a new link. Only one circumstance prevent their en- 
richment - seizure of this link by another node. In such 
kinetics of distribution of links, there are no any finite ra- 
dius of "interaction" and there are no principal obstacles 
for capture of great fraction of links by some node. 



r / "exp[-/3Ae 3 F 2 (l,l,l + A;2,2;0] , (40) 

where C = 0.5772 ... is the Euler's constant and ip( ) is 
the ■i/'-function. The transcendental equation for /3 may 
be written if one substitutes Eq. ( (40| ) into the right side 
ofEq. @: 



/r 1 



e _ /3 [C+^(l-A)] x 



o C^(i-0 A 



exp [-f3A( 3 F 2 (1, 1, 1 + A; 2, 2; C)] • (41) 



(An equivalent equation one may obtain substituting the 
solution, Eq. (^0|), into the relation, f Q d(n(() = 2). The 
7 exponent may be obtained using the universal relation, 
Eq. (H). 

The solution of Eq. ( [4l| ) exists in the range, — oo < 
A < 1. One may also find simple expressions for /3(A) 
and 7(A) at A — > 0: 



/3 = - — (1 — In 2) A , 7 = 3 + 4(1 - In 2) A . 



(42) 



where the numerical values of the coefficients are 1 — 
In 2 = 0.3069... and 4(1 -In 2) = 1.2274.... We used the 
relation 3 F 2 (1, 1, 1; 2,2; C) = Li 2 (C)/C = (EZi C k /k 2 )/( 
while deriving Eq. (42). Here, Li 2 ( ) is the polylogarithm 
function of order 2. 

In the limit of A — > 1, using the relation 
3 F 2 (1, 1, 2; 2, 2; C) = - In(l - ()/(, we find 



/3 = Cl (l-A) , 7 



1 1 
ci" 1 — A 



(43) 



1.2400. 



the constant c\ 



Here, c\ = 0.8065 . . . , cj~ 
is the solution of the equation 1 + 1/ci = exp(ci). The 
dependences, /3(A) and 7(A), are shown in Figs, g and^. 
For comparison, data of simulations pl| are also plotted. 
One sees, that (3 — > 1 and 7 — > 2 in the limit A — » —00. 
Therefore, the whole range of the variation of j3 is (0, 1) 
and of 7 is (2, 00). One should note that, in this section, 
we studied the case of m = 1 but the results, the scaling 
exponents are, of course, independent of m. 



VI. CAPTURE OF LINKS BY "STRONG" NODES 

Now we can consider effects of s-depending strength 
of nodes, G = G(s). We proceed with our tactics and 
demonstrate these effects using minimal models. 

It is worth to start this section from the following no- 
tion. All networks that we study in the present paper 
have a general feature - each their node has a chance to 



A. Varying strength of nodes 

Let us start from the case of a power-law dependence 
of the strength, G(s) cx s~ A . At first sight, this problem 
seems to be very close to the one that was considered in 
the previous section. Nevertheless, repetition of calcula- 
tions of Sec. [y] for the present model leads immediately 
to the following expression for the degree, 



-A'l 



(44) 



where B is a constant, £ = s/t, k = k(£) (compare with 
Eq. (38)). This expression provides quite different be- 



havior of k(s,t) than that one that we considered in the 
previous section. 

If A' < 0, i.e., new nodes are "stronger" than old ones, 
we get n(s/t — > 0) — > 0, and network is exponential. If 
A' = 0, we get an ordinary scale-free network. If A' > 0, 
then K,(s/t) is extremely divergent function at £ = 0. This 
peculiarity is not integrated. It indicates that several old- 
est nodes capture a finite part of nodes. Our theory is 
not applicable in this case directly, since one can not 
solve a self-consistency equation and find f3. We do not 
consider this situation in detail here. Such behavior may 
be compared with that one found for G = G(k) = fc^ -1 , 
y > 1, where most of links turn to be attached to the 
oldest node [|l6|,E5). The "phase diagram" of the model 
is shown in Fig.Ta 



B. Local strong nodes 

Fluctuating strength of nodes, G(s) is introduced in 
the resent papers of Bianconi and Barabasi, Ref. p4] , p2"|| . 
Using analogies with other cooperative problems of sta- 
tistical mechanics, they indicate a new phenomenon aris- 
ing in this situation. We, on the contrary, demonstrate 
the essence of evolution of such systems using the most 
simple (although quite general) example, without imple- 
menting any analogy. 

Let us start from the model with directed links intro- 



duced in Sec. IV B. To simplify the formulas, we assume 
that A — (one may see that this does not reduce ex- 
tend of generality of the model which produces scaling 
exponents in the wide ranges of values, 2 < 7 < 00 and 
< (3 < 1). Let us recall the model (see Fig. |l|, a). 
Each increment of time, a new node with n attached to 
it links is added to the network. These n links are di- 
rected to the new node. Simultaneously, m extra links 
are distributed preferentially amo ng no des. The rule of 
preference is the same as in Sec. IV B, i.e., probability 
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that a link is directed to a node s is proportional to its 
in-degree, q s but with one exception - one node, s, is 
"stronger" than others, that is probability that this node 
capture a preferentially distributed link is higher. It has 
an additional weight factor, g > 1, and proportional to 
gqs- This means that G s = 1 + (g — l)d s j. 

Let us study behavior of the network at t ^> s. The 
equations for average in-degree are 



gm 



dq- a (t) 

dt 

dqjs, t) 
dt 



9Q&) 



(.9 - + Jo dsq(s,t) 



q(s,t) 



(.9- l)q s (t)+ f dsq(s,t) 



> <?s(* = s) = Qi, 
, q{t,t)=n. (45) 



Here, + J dsq(s,t) — (m + n)t + O (1) . One can see 
that this equality alone is not sufficient to determine the 
denumerators in Eq. (jig). 

Two different situations are possible. In the first one, 
q$(t) grows slower than t, and, at long times, the denom- 
inators are equal (m + n)t. Then, second line of Eq. j45| ) 
is similar to Eq. ([22]), and we get again the exponents, 
P = raj (m + n) = /3q and 7 = 2 + n/m = 70 = 1 + l//3o> 
were < /3q < 1, 2 < 70 < 00. (It is convenient to 
write equations using not the parameters, m and n, but 
70 or Pq, that is the exponents of the network in which 
all nodes have equal strength, g = 1.) The first line of 
Eq. (Sq), in this case, looks as 



dq~ s (t) _ gm gq s (t) 



dt 



m + n t 



(46) 



Hence, at long times, q s (t) = const(qi)t gm ^ m+n \ and 
we see that the in-degree of the strong node does grow 
slower than t only for 



g < g c = 1 + — = 70 - 1 = P 1 > 1 
m 



(47) 



so we obtain the natural threshold value. 

In the other situation, g > g c , at long times, we have 
the only possibility, qs(t) = dt, d is some constant, since 
more rapid growth of qg{t) is impossible, d < m + n. 
This means that, for g > g c , finite part of all preferen- 
tially distributed links is captured by the strong node (in 
p2[ this situation is called condensation). We see that 
a single strong node may produce a macroscopic effect, 
so g c is a giant capture threshold. In this case, Eq. J45| ) 
takes the form, 



gm 



?.-(*) 



dUt) = 

dt (g-l)d + m + n t ' 

dq(s,t) m q(s,t) 



(g — l)d + m + n t 



(48) 



Note that the coefficient in the first line is always larger 
than the one in the second line, since g c > 1. From the 
first line of Eq. (Est), we get the condition 



(g — l)d + m + n 



(49) 



so, for g > g c , the following part of all links of the net- 
work is captured by the strongest node: 



1 9- 



9c .9-1 



(50) 



We have to emphasize that d is independent on initial 
conditions! (Recall that we consider the long time limit.) 
This giant capture leads to change of exponents. Using 
the condition Eq. ( f49| ) , we get immediately the following 
expressions for them, 



P=-<Po- 
9 



7 = 1 + g > 70 . 



(51) 



The fraction of all links captured by the strongest node 
and the P and 7 exponents vs g are shown in Fig. |^. Note 
that growth of g increases the value of the 7 exponent. If 
World is captured by Bill Geits or tzar, the distribution 
of wealth becomes more fair! One should note that the 
strong node does not take links away from other nodes 
but only intercepts them. The closer 70 to 2, the smaller g 
is necessary to exceed the threshold. Above the thresh- 
old, values of the exponents are determined only by g. 
Nevertheless, the expression for dj (m + n) contains 70 or 

For g > g c , in the giant capture regime, the strongest 
node determines the evolution of the network. With 
increasing time, a gap between the in-degree of the 
strongest node and maximal in-degree of all others grows, 
see Fig. ^. A small peak at the end of the continuous 
part of t he dist ribution is a trace of initial conditions, 
see Sec. [VIII B . Note that the network remains scale- 
free even above the threshold, i.e., for g > g c , although 
7 grows with growing g. 

The same result may be obtained for several nodes, §i 
with different strengths gi, where Si -C t. Here, <?; > 1 
(<7i < 1 do not produce any visible effect). In this case, 
the strongest node, Sj, again captures a finite part of 
nodes if gj > g c . Note that the time which the strongest 
node needs to seize a finite fraction of links may be very 
long if it is only a bit stronger than the previous strongest 
one. Note also that a death of the strongest guy produces 
a dramatic effect - after some time, the distribution be- 
comes less "fair". Deleting of a single strong node may 
also destroy entire network. 

If N strong nodes of equal strength, g, are present, one 
can find that, again, the critical exponents are described 
by Eq. ( pl| ) and the part of links captured by all these 
strong nodes together is given by by Eq. (p0[) ■ 
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C. Mixture of "weak" and "strong" nodes 

How one may "smear" the giant capture threshold de- 
scribed in the previous section? Again, let us use the sim- 
plest model in which nodes have two values of strength, 
1 and g > 1, but now the nodes with different strengths 
are distributed homogeneously. Probability that a node 
has the strength equal 1 is 1 — p, and, with probability p, 
a node has the strength g. All other conditions are the 
same. We can introduce average in-degrees, <Zi(i) and 
q g (t) for these components. Then, the evolution of the 
average in-degrees is determined by the following equa- 
tions, 



9Q q (s,t) 



(l-p)f dsq 1 (s,t)+gpf ds q g {s,t) 



2i(s,i) 



(52) 



dq g (s,t) 
dt 

dq t (s,t) 

5 = 1 1 

(1-P)/o ds<h(M) +9Pf dsq g (s,t) 

(compare with Eq. (ESI)). As usual, we get 



(1 — p) / dsq^s, t)+p / ds q g (s,t) — (m + n)t . (53) 
J o Jo 

If one introduce the natural notation, m/[(l — 
P)Jo dsq^s^) + gpj*dsq g (s,t)} = f3 u Eq. @ takes 
the following form, 



dt 



= Pi9- 



dq~i(s,t) qi(s,t) 
= Pi- 



Ot 



t 



(54) 



Inserting the solutions of these equations, q g (s,t) = 
n{s/t)-^ a and q x {s,t) = n(s/t)- pi into Eq. (§f), we 
obtain the equation for /3i, 

1 -p p 



1 



1-/3! 1 - fog l-p 
The solution of Eq. (|55|) is of the form, 



(55) 



0i = ^-{[l-p+pPo} + \p+(l-p)0o}9- 
2.9 

^{[i-p + pPo] + \p+(i- p)Po]g} 2 - Wog } (56) 

(we chose this root of the square equation to obtain the 
proper equality, 0i(g = 1) = /3o). 

The fraction of links captured by the strong component 
is given by the relation, 



d 



P 



i-A) 



m + n 1 — Pig m + n 1 — j3\ 



(57) 



with /3i taken from Eq. (^6|). Exponents (3\ and 71 = 
1 + l//3i describes the distribution of nodes of the weak 



component containing the 1 — d/(m + n) part of nodes. 
Exponents /3 g = Pig and j g — 1 + l//3 g describes the dis- 
tribution of nodes of the strong component. The depen- 
dences of these characteristics on g are shown schemati- 
cally in Fig. |[ One may see that, for p — > these curves 
tend to the dependences obtained for a single strong node 
in the previous section. For p > 0, a long tail of the dis- 
tribution with the exponent 7 ff < 70 < 71 is determined 
by strong nodes. At p — > 0, it is transformed into the 
(^-function obtained in the previous section. A particular 
value of p determines the smearing of the giant capture 
threshold. 



VII. DISTRIBUTIONS OF LINKS BETWEEN 
PAIRS OF NODES 

Let us discuss another characteristic describing a struc- 
ture of networks - distribution of links between pairs of 
nodes - an average matrix element of an adjacency ma- 
trix. (An element b% j of an adjacency matrix equals 1 if 
there is a link connecting the sites i and j and it equals 
is the link is absent.) 

There is a very important particular case, when this 
characteristic may be may be easily calculated. In this 
case, new links appear only between new node and old 
ones but never - between old nodes of a network (see 
Fig. [l], b). For instance, networks of scientific citations 
belongs to class of such networks. 



Here, we study the simplest model of Sec. IV B with 
n = 0. For brevity, we consider in this section net works 
with undirected links, so, into the relations of Sec. IV B 
we may substitute q(s, t) = k(s, t) — m. That is, at each 
time step, we add a new node with m attached undi- 
rected links. Other ends of these links are distributed 
preferentially between old nodes. 

Let us introduce the average number of links between 
nodes s and s' at time i, b(s,s / ,t). For a network that 
we consider, 



k(s,t) = J dub(u,s,t) + J dub(s,u,t) — 

m+ dub(s,u,t) . (58) 

J s 

No new links between old nodes arise, so one can see that 

dk(s t) 

b(s, t, t'>t) = b(s, t, t) , — = b(s, t, t' >t). (59) 

Using scaling representation, k(s,t) — n(s/t) and 
b(s, s' , t) = t~ 1 B(s/t, s'/t), one can write 

K»,s', t '>0 = >(ii) = i B (j4 = i) (60) 

and 
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d{s/t) 



«(?) 



d{s/t) 



t 1 1 



(61) 



Therefore, we get the following relations for #(£,£') and 



s(e,i) = -e«'(0 



(62) 



where k'(ie) = dn(x)/dx. 

In our particular problem, k(1) = m, so, from Eq. 
([58|), it follows that 



«(0 



J d?'BW) = m + J jB(C>1) 



From relations of Sec. [V B| , we obtain 



«(0 = 

so we get the result, 



1 



B(e,0 = m(i-j8)r /, r / '- 1 



(63) 



(64) 



(65) 



This characteristic was obtained explicitly for a model 
with (3 = 1/2 p6|| . From this exact solution, the par- 
ticular case, /3 = 1/2, of Eq. ( |65| ) follows immediately. 
Note that £?(£,£') does not proportional to the product, 

It is interesting to compare these expressions with the 
corresponding results for exponential networks (we put 
off the discussion of applicability of the continuous ap- 
proach to non-scale-free networks until Sec. [X Q ). Set- 
ting m — n r and n = A 
one gets 



in the equations of Sec. [IV C 
k(0 = m(l - InO , B(£,0 



m 



(66) 



Eq. (J6e|) follows from Eqs. fl64|) and fl65j) at /3 0. 

We have obtained the density of linkage between pairs 
of nodes with fixed times of birth. From this charac- 
teristic (in the continuous approximation), we can find 
average number of connections, T>{k, k'), of parent nodes 
with degree, k, and child nodes with degree, k'. In the 
continuous approximation, k > k' . 

Let us consider only the scaling region and will not 
write coefficients in the following formulas of the present 
section, k = so, for a scale-free network, P{k) = 

-{dn/d^ 1 ~ fc~ 1//3_1 . Analogously, 

V(k,k') = 



k-^k'- 2 . 



One sees that T>(k, k') is decreasing function of both k 
and k' . Therefore, if k (parent node) is fixed, the most 
probable linking is with a child node with the smallest 
possible k', i.e., with the node with k' ~ 1. 

If k' (child node) is fixed, the most probable linking 
is with a parent node with the smallest possible k, i.e., 
with the node with k = k! . This answer may be com- 
pared with the corresponding exact result of Krapivsky 
and Redner |35f| - the maximum of the probability that 
a node with degree k (parent) and a node with degree k' 
(child) are connected, when k' is fixed, is at k/k' — 0.372. 
Therefore, we see, that the continuous approach is good 
also in this situation. 

Eq. (^) looks asymmetrically and cannot be factor- 
ized. The reason is the obvious absence of the symmetry 
relatively reversal of time - a quite natural asymmetry 
between parents and childs. 



VIII. ACCELERATING GROWTH 



Above we studied only linear network growth, i.e., the 
total number of links of a network divided by the total 
number of its nodes was constant during the evolution. 
This is only a particular case of network growth. The 
total number of links may be a nonlinear function of the 
total number of nodes. Keeping in mind the most in- 
triguing applications, communications networks, we con- 
centrate now on accelerating growth |33| where this func- 
tion grows more rapidly than a linear one. In particular, 
we will see that a power-law dependence of input flow of 
links produces scale-free networks. 



A. Scaling relations 



In the present section, we consider scale-free networks 
with power-law dependent on t input flow of links. We 
will see that, in a large network size limit, such nonlinear 
growth may produce the following non-stationary distri- 
butions, 



and average indegree (or degree), 



q(s, t) cx t 6 



(68) 



(69) 



Let us discuss what is the result of such non-stationary 
behavior. 



(67) 



One can easily repeat the derivation of Sec. Ill and 
obtain the scaling form of the connectivity distribution 
of individual nodes: 
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) / 




:d1 



(70) 



We again use the definition of P(k,t), Eq. (m), and Eq. 



dxt^x 13 }{qt~ S x P ) cx l*/P q -l-W cx tV 7 . (71) 



Hence, we get the universal scaling relation for the expo- 
nents z, S, and f3, 



z = 5/p 



(72) 



and the old relation, (|l2]), for the exponents 7 and f3. 

Let us derive this result once again applying the con- 
tinuous approach. Using Eqs. (pd]), (|68|), and (|69|), we 
obtain 

f g -7 oc -~[q = t s +P S -P] cx ^g-i-V/J . (73) 



From Eq. @, the relations, Eqs. Q and (|7|) follows 
immediately. 

One can write the relation, Eq. (|70| ) using only the 7 
exponent: 



p(q,s,t) 



,1/(7-1) 



4(1+^/(7-1)^ i 9 4(1+^/(7-1) 



,1/(7-1) 



Eq. ( |74| ) is direct generalization of Eq. ( |14[) obtained for 
a linearly growing networks. 



B. Cut-off of connectivity distributions 

Eq. (|68| ) is valid only in the limit of large network 
size (long times). Let us discuss briefly arising finite-size 
effects. 



Account for the relation Eq. (|69|) of Sec. VIII A to- 
gether with the rapid decrease of p(q, s, t) at large q pro- 
duces a cut-off of power law distributions at the charac- 
teristic value, 



lent 



t (3+S = t P(l+z) = t (l+*)/( 7 -l) 



(75) 



Here, we used the fact that q(s, t) is largest for the oldest 
node and the scaling relations between the exponents, 
Eqs. ( |l2| ) and ([72]). It was shown that a trace of initial 
conditions at q ~ q cut may be visible in degree- and in- 
degree distribution measured for any network sizes p6[ . 
Such a cut-off (and a trace of initial conditions) inputs 
restrictions for observations of power-law distributions 
since there are few huge networks in Nature. 

Let us obtain a general form of P(q, t) for scale-free 
networks in scaling regime. Using the known scaling form 
of p(q, s, t), Eq. (|74|), we can write, 



P(q,t) 



dsp(q,s,t) 



t z q-t / dww^^fiw). 

y g t-(i+»)/(-r-i) 



(76) 



Passing to the scaling limit, q — * 00, t — > 00, and 
q,£-(i+z)/(7-i) j g f} xeclj we can replace the upper limit of 
the integral in Eq. (|7^) by 00. Then, we get immediately 
the following scaling form, 

P(g,i) = t z q-~>F(qt-( 1+z ^) = t*g-T.F(gi-(l+*)/(T-l) ) , 

(77) 

where F(w) is a scaling function. In the case of linearly 
growing network, the exponent z is zero, so P(q,t) = 
q-^F{qt- 1 /P) = g-TF(gt- 1 /(7-i)). This relation was 
obtained for an exactly solvable model p™ . 



C. Scaling exponents 



VIII A and VIII B 



Using relations obtained in Sees 
one can obtain general results for a nonlinear, accelerat- 
ing growth of networks. 

We start from the most general considerations. In 
the scale-free networks, a wide range of the connec- 
(74) tivity distribution function is of the power-law form, 
P(q) oc t z q~~<. It will be clear from the following that, 
to keep a network in the class of free-scale nets, the flow 
of new links has to be a power function of the number of 
nodes of a network, i.e., be proportional to t a . 

First, let us assume that the exponent of the distribu- 
tion is less than two. The reasonable range is 1 < 7 < 2. 
To produce the restricted average connectivity (that is 
proportional to t a ), the distribution has to have a cut 
i(i+*)/(7- ' 



< 



(see Sec. 

q c . The restriction 



VfflB| ). 



off at large q, q c 
P(q) oc t z q-i for q t - t x < q 
from below is necessary to garantee convergence of the 
integral, J °° dq P(q,t) — 1. From this, we get immedi- 
ately x = 2/(7— 1). (Of course, this relation is valid also 
for 7 > 2.) 

The average connectivity, q is of the order t a+1 /t, then 
t a - / dqqt z q~~< - f-i+(i+*)/(7-i) (the value 

of the integral is determined by its upper limit). There- 
fore, (1 + z)/(-f — 1) = a + 1, so the cut-off of the distri- 
bution is of the order of the total number of links in the 
network. This is the maximal number of the problem, 
hence, in fact, the cut-off is absent. The expression for 
the 7 exponent, 7 = 1 + (1 + z) /(l + a), follows from the 
last relation. 

Note that a is "external" exponent which governs the 
growth process. We have demonstrated that it is enough 
to know a and only one exponent of 7, /?, z, 5, x for finding 
all the others. 

Note also that we have to set z < a to keep the expo- 
nent 7 below 2 as it was assumed. Also, one sees that 
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the lower boundary for 7, 1 + 1/(1 + a), is approached for 
the stationary distribution, z = 0. In this case, the form 
of the distribution is completely fixed by the accelerating 
growth, the exponent 7 depends only on a. 

The other possibility is 7 > 2. In this case, the in- 
tegral for the average connectivity is determined by its 
lower limit, t a ~ f t */^-i) dq qt z q~f ~ t*-*(7-2)/(7-i). 
Hence, 7 = 1 + z/a, and z > a. (Of course, this rela- 
tion is not valid for a = 0.) Thus, we have described the 
possible forms of the connectivity distribution. 

Let us demonstrate how these distributions may arise 
in the nonlinearly growing networks with the preferential 
linking. We int roduc e the simplest generalizations of the 
model of Sec. IV B to the case of the increasing input 
flow of links, cot a , a > 0. 

First, let us consider the case of constant initial attrac- 
tivity, A — const. The equation for q{s,t) is 



dt 



cot a 



q{s,t)+A 



J du [q(u,t) + A] 



(78) 



g(0,0) = 0, q(t,t) — n. One may check that 
L duq(u,t) = nt + cot a+1 /(a + 1). Inserting this re- 
lation into Eq. (|78[ ) and solving the resulting equation, 
one gets 



q{s,t)+A 
n + A 



l + (n + A)(l + a)t- a /c 
1 + (n + A)(l + a)s- a /c a 



1+1/ a 



S\ -(<*+!) 



(?) 



In the interval [(n + A)(l + a)/c ] 1/Q < s < t, 



q(s,t) = (n + A) (|) 



(79) 



(80) 



Thus, the exponent f3, q(s,t) cx s 13 , equals 1 + a and is 
larger than 1. The dependence q(s) becomes constant, 



at s < \{n + A){\ + a)/co] 1 ^ Q . One may compare the 
result, Eq. (|8|), with the total number of links in the 
network, N(t)& c t a+1 /(l + a). 

From Eq. (|8l|), we see immediately that the exponent 
/3 = l + a, so 7=1 + 1/a, d = z = 0. One may cal- 
culate the connectivity distribution using Eq. (pi). The 
resulting distribution, in the region 1 <C q/(n + A) <C 
{ c /[(n + A)(l + a)} } 1+1/Q is of the form, 



1 + a 



(82) 



Thus, we obtain the stationary connectivity distribution 
with the 7 exponent less than 2 that belongs to one of 
the types described above. 

To demonstrate the other possibility, 7 > 2, we con- 
sider below the model with a different rule of the distri- 
bution of new links. Let the additional attractivity be 



time-dependent, and new links be distributed between 
nodes with probability proportional to q + Bcot a / (1 + a), 
where B is positive constant. cot a /(l + a) is the average 
connectivity of the network at time t. 

Repeating the previous calculations, one gets the equa- 
tion 

dq{s,t) q{s,t)+Bc t a /(l + a) 

— Co t 



dt 



nt + Bc £"+7(1 + a) + c t a+1 /(a + 1) 



g(0.0) = 0, q(t,t) = n. At long times, one obtains 

dq(s,t) _ 1 + a q{s,t) + Bc t a /(l + a) 
dt ~ 1 + B t 

The solution of Eq. (pj) is 



(83) 



(84) 



q(s,t) = 



Bc s a 



1-Ba 



S x _(l+a)/(l+B) Bc t° 



(?) 



1-Ba 



(85) 



If B = 0, we obtain the previous result, (3 = 1 + a. For 
s q »!i(1-5q)/(Bc ), 



q(s,t) 



Bc t a ( /s\«-(i+°0/(i+s) 



1-Ba 



(?) 



(86) 



Therefore, the scaling exponents of the growing network 
are/3= (1 + a)/(l + B) - a = (1 - Ba)/(1 + B), 7 = 1 + 
1/(3 = l + [{l+a)/{l + B)-a]- 1 = 2 + B(l + a)/{l-Ba), 
S — a, and z = a(l + B)/(l — Ba). The connectivity dis- 
tribution differs sharply from the distribution obtained 
for the previous model. It is nonstationary and is of 
the form P(q,t) ~ t ^ 1 + B y^- Ba )q-^ 1+B )^ 1 - Ba ^ for 
g> t a . In this case, (3 <1 and 7 > 2 for any positive a 
and B. The scaling regime is realized if Ba < 1. 

Note that, in both considered cases, one cannot set 
a = directly in the obtained expression for the scaling 
exponents. 



IX. DEVELOPING AND DECAYING 
NETWORKS 



Now we can study evolution of network accompanied 
by reconstruction of its old part ]32[ ]. This, e.g., may 
include permanent deletion of old links or nodes. Note 
that the processes of addition and deletion of links may 
be considered in a unified way, so we study them to- 
gether. For brevity, in the present section, we consider 
undirected links. 
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A. Developing networks 

Let us introduce two channels of app earan ce of new 
links. The first one was studied in Sec. IV B| . We con- 
sider undirected links, so, in our old formulas, we have 
to substitute q(s,t) = k(s,t) — m and put k(t,t) = m. 
n = 0. Instead of the additional attractiveness A we 
use here the constant A n = A — m > —m. The second 
channel is the following. Each time a new node is ap- 
peared, c additional links arise between old unconnected 
nodes i and j with probability proportional to the prod- 
uct (ki + A a )(kj + A a ), A Q > —to. Note that A n and A Q 
may be not equal, so, here, we h ave a mixture of differ- 
ent preferences like in Sec. IV C| . Then the equation for 
k(s,t) is of the form, 



dk(s,t) 



k(s, t) + A„ 



2c- 



k(s,t)+A D 



9t J Q du[k(u,t)+A n ] ' ' f*du[k(u,t)+A ] 

(87) 

(see [32]]). J Q ds k(s, t) = 2(m + c)t, so we obtain imme- 
diately the scaling exponent (3, 







2c 



2m + 2c + A n 2m + 2c + A Q 



(88) 



and using the relation Eq. (|T^), the exponent 7. We do 
not consider these relations in detail. 



B. Decaying networks 

Let us introduce possibility of deleting of links. The 
second channel now is the following. Each increment of 
time — c random links are removed from the network. 
Note that we we save the same definition of the c as in 



Sec. 1XA to use similar relations. 



Previously, a process of instant random damage was 
considered [ p||l5| , p4"| ,?]. Such type of damage can not 
change a value of 7 neither for random removal links nor 
for removal nodes. Here, we consider a quite different 
situation, permanent random damage - components of a 
network are removed permanently during its growth. In 
this case, difference between random removal of nodes 
and links is striking. 

We have shown that Eqs. ( |87| ) and ( |88| ) also describe 
the case of decaying networks with permanent deletion 
of links if one takes c to be negative in them and assumes 
that A = §1. Then, it follows from Eq. (||) that 



1 



+ 



c/m 

2 + 2c/m + A n /m ' 1 + c/m 

1 + 7oc/to + 2(c/m) 2 
(1 + c/to)(7o-1 + 2c/to) 



(89) 



7 = 2 



1 + 7oc/m + 2(c/m) 2 
70 - 2 + c/m 



(90) 



where 70 = 7(0 = 0) = 3 + A n /m. 

The resulting phase diagram, —c/m vs 70 is surpris- 
ingly rich, see Fig. One sees that the effect of deletion 
of links during the network growth may be quite different 
for different values of the 70 exponent of the correspond- 
ing network without such deletion. For 70 < 1 + v2, the 
removal of links decreases 7. If 1 + \[2 < 70 < 2\f2, 
deleting of links first increases 7 but afterwards 7 de- 
creases. For 70 > 2\f2, 7 grows monotonously with 
increasing —c/m until it becomes infinite on the line, 
70 = — 2(c/m) — 1/(c/to). We have checked our results 
by simulation in the particular case of 70 = 3 [ p2| . 

One should note that permanent deletion of links pro- 
duces a region between the lines, — c/m = (70 — l)/2 and 
— c/m = 7o — 2, in which 7 varies between 1 and 2. It is 
the first time we find these distributions, with 1 < 7 < 2, 
in such problems. In dashed regions of Fig. ||, station- 
ary distributions are absent. Note that, for large enough 
—c/m, the network may decay for a set of uncoupled 
clusters. 

Let us consider now additional permanent deletion of 
random nodes. Ab initio, one may think that this factor 
does not change 7. Nevertheless, as we shall see, this case 
is of a special interest. 

At each increment of time, one new node is added, 
and a randomly chosen node is deleted with probability 
c < 1. We use the preferential linking introduced in Sec. 
IV B, although a more general model may be also con- 



sidered. It is natural to consider in-degree, q(s,t), here, 
since results are more general in this case. 

Let us introduce the probability that a nodes s is 
present, A/"(s, t). In the continuous approach it has mean- 
ing of "density" of surviving nodes at time t , J\f(s, t) . One 
may introduce "density of in-degree" . In the continuous 
approach, it looks as 

d{s,t) =J7(s,t)q(s,t) . (91) 

Our main equation for J\F(s, t) is of the form, 



dW(s,t) Ms,t) 
- — c - 



dt 



J duM(u,t) 



(92) 



where 77(0,0) = and 77{t,t) = 1. From Eq. (||) one 
gets immediately the obvious relation, L duU(u,t) — 
(1 - c)t. The solution of Eq. (||) is 

A7(M) = (-J ■ (93) 



and 



In the present case, the equation for the in-degree is 

(94) 



dq(s, t) q(s, t) + A 

— m 



dt 



Jo duJ7(u,t)[q{s,t) + A] 
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where the integral on the right side is equal to the sum of 
the total in-degree of the network at time t and the prod- 
uct of the additional attractiveness and the number of 
sirvived nodes (compare with Eq. (|l9| ) for c = 0). Here, 
5(0, 0) = and q(t, t) = n so d(0, 0) = and d(t, t) = n. 

Multiplying Eq. J94|) by 77(s,t) and applying f* ds to 

both sides of the resulting equation, we get J* ds d(s, t) = 
(to + n)(l — c)t. Substituting this relation into Eq. (94) 



we get immediately q(s, t) oc (s/t) 13 , where the exponent 



00 



(1 - c)(m + n + A) 



1 



(95) 



Here, /3 = (3(c = 0), 70 = 7(0 = 0). Using Eq. 
(|93|), we obtain for the density of in-degree, d(s,t) oc 
(s7*r (/3+c)/(1 ~ c) - Note that, d(s -> 0, i) -> 00 for c < /3 , 
and d(s — > 0, t) — > for c > /?o- 

The expression for the degree distribution looks as 



P(k,t) 



f*dsN(s,i)p(q,s,t) _ 
J*dsAf(s,t) 

J*ds7T(s,t)6(g-q(a,t)) 
J*dsjf(s,t) 



(96) 



III 



Therefore, repearting the derivation of Sec. 
we obtain the following distribution, P{k) oc 

fc-(l/fl[c/(l-c)]fc-l-l//S fx fc-l-l/[/3(l-c)] ; so 7 = 1 + 

l/[/3(l - c)] = 1 + 1/A> = 7o- Thus, the distribution 
is of the same form as without the permanent random 
deleting of nodes. Note that there is a great difference 
from the case of the permanent random deletion of links 
considered above. 

Here, we find for the first time a violation of the scaling 
relation, Eq. (||). The reason is effective renormaliza- 
tion of the s variable due to the removal of nodes. 

Repeating scaling considerations of previous sections, 

we get for this case the following forms of p(q, s, t 

and P(q), p(q,s,t) = {s/tf f[q{s/tf) and P(q) = 
q -{W/Wi-c)]} F ( q / t Py 

We conclude this section by the statement that perma- 
nent deleting of a part of nodes with the largest values of 
degree, that is an analogy of intentional attacks |^3|,^5| , 
destroys scaling behavior of a network. One may easily 
check this statement using the continuous approach. 



X. APPLICABILITY OF THE CONTINUOUS 
APPROACH 



In the present section we discuss quality of the used 
continuous approach. For this, we compare known exact 
results with results obtained in frames of the continuous 
approximation. 



A. Linear preference 

In Sec. |l[ we have written out already the answer of 
the continuous approach for the Barabasi- Albert's model 
at m = 1. For a more general case when m is any posi- 
tive integer number, and k = q + to, in the frames of the 
continuous approach, one gets, P(q) = 2m 2 /(q + to) 3 
|l8| . One may compare this expression with the ex- 
act result obtained without passing to continuous limit: 
P(q) = 2to(to + \)/[{q + m){q + to + l)(q + to + 2)] @. 
Exponents are the same but the factors are different. 

One may ask, why the approach is so good. The reason 
is rapid decrease of p(q, s, t) at large q. Because this, the 
results of the continuous approximation obtained with 
the S- function ansatz, Eq. (^), are reasonable. Indeed, 
in the limit of large s and t and for fixed s/t, we obtained 
the following scaling exponentially decreasing expression, 
p{q,a,t) = \q(s/m A -i e^p\-q(s/m/T(A) @. The 
model of Sec. IV B with n — was considered. Here 



r( ) is the gamma- function. 



B. Nonlinear preference 

The continuous approach may be applied to nonlinear 
preference. In this case, calculations similar to that ones 
of Sec. may be made. 

Let us consider the simplest case, generalizing the 
Barabasi- Albert's model, k(t,t) = l,m = 1. Then, the 
main equations are, 



dk(s,t) 
dt 



f P [k(s,t)] 
J Q du f p {k(u,t)\ 



dsk(s,t) = 2t. (97) 



Here, f p (k) is a preference function. Let us search for the 
solution of Eq. ( |97j ) in the scaling form, k(s, t) = K(s/t). 
Then, 



fp["(0] 



91n ^ lodCfMO] 



k(1) = 1 



dC«(C) 



(98) 



Let us start from the case of nonlinear preference that 
produces scale-free networks. Let the probability for dis- 
tribution of new links be proportional to the preference 
function, f p (k), that is linear asymptotically for large k, 
f p (k — > 00) — > ck, where c is a constant, (see |l^ , ^5|] ). 

The integral J* d( / p [k(C)] is a constant of the problem. 
In the scaling region of large k(s,t) and k, the equation 
takes the form, 



-dln«(0/dln£ = c|y <Kf p [K(Q]j =(3. (99) 

Eq. ([39]) demonstrates that the scaling is present, and 
the network is scale-free in this case. 
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Therefore, to find the scaling exponent f3, we have to 
solve the equation, 



(100) 



and after inserting its solution, «(£) = F _1 [F(1) + 
(/3/c)ln£] (F(k) = J dnf^n), F" 1 is an inversed func- 
tion to F), into 



/3~ 1 =c- 1 d(f p [K(()} or equivalents 2 = / d( k(() 
Jo Jo 

(101) 

find the solution f3 of any of these transcendental equa- 
tions. 

We will not consider examples of application of these 
relations but describe briefly the case of a power law pref- 
erence function just to test the continuous approach us- 
ing the known non-scale-free network. Indeed, we have 
checked the quality of the continuous approach for scale- 
free networks. Now it is naturally to check it for other 
networks. 

Let the preference function be f p (k) = k~ v . If one sets 



J d(^K y (() = ji = const , 



(102) 



then the equation for is 



- dln ^' = /x- 1 exp[-(l-2/)ln«]. (103) 



Its solution is 



A 1 



(104) 



The constant \x can be obtained from the transcendental 
equations 



which gives /i(y). Near y = 1, /i = 2y, near y = 0, 
A* = 1 + 0.5963y, where 0.5963 = eEi(-l). 

Inserting the solution, Eq. (104), into the expression 
for P(k) in the continuous approach, Eq. ( |TT| ) we obtain 
the connectivity distribution, 



P(k) = 



y) fc-y 



exp 



-k 1 - y 



(108) 



These results are close to the exact ones [^6|35J]. The 
values of the powers are the same, although the coeffi- 
cients, /j differs slightly. E.g., in Refs. |l(||35| Near y = 1, 
H 2 - 2.407(1 - y), and, near y = 0, /i S 1 + 0.5078y. 



C. Random linking 



Finally, in the frames of the continuous approach, let 



us consider the model of Sec. [VA producing non-scale- 
free, exponential degree distribution. This network is a 
continuous case of the network considered in Sec. 
n = 0, m = 0, n r = 1, A = 0. The solution of Eq. (| 
this particular case is 



IVC, 



k(s,t) = 1 - ln(s/t) 



(109) 



(we used the boundary condition k(t,t) = 1). Then, us- 
ing Eq. (pi]), one gets the connectivity distribution, 



P(k) = - 



1 <9[texp(l - fc)] 



dk 



exp(l — fc) . 



(110) 



IV A 



This can be compared with the exact form of Sec 
P(fc) = exp(-fcln2). 

Therefore, the continuous approach easily produces 
reasonable answers even for non-scale-free networks. The 
reason is again rapid decrease of p( fc, s, t) at large fc for 
these networks, see Eq. ( |l8| ) of Sec. |1V A 



CONCLUSIONS 



2 = / d£ I 1 - -ln£ 



1/(1-1/) 



or equivalently 



o 



y/a-y) 



(105) 



(106) 



The final transcendental equation may be written in 
the form: 



2 1^-1 e*)=r 1+ ' 



1-1/ 



i - y i - v 



(107) 



We have analysed scale-invariant properties of scale- 
free networks with growth governed by a mechanism of 
preferential linking. Degree distributions of such net- 
works are of a simple scaling form. We have shown that 
the arising scaling exponents are coupled by universal 
scaling relations. Nevertheless, we present an important 
particular case in which these simple relations arc violate. 

One of discussed questions has been about types of 
preferential linking producing scale-free networks. One 
can see that scale-free networks are produced by wide 
variety of linking. In particular, it is enough to add an 
admixture of linear preference linking to random linking 
to get a scale-free network. Interplay of different factors 
such as deletion of links of a network during its growth 
may change dramatically its degree distribution and even 
remove it from the class of scale-free networks. Hence, 
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we have shown, how one can increase or decrease critical 
exponents of a network. 

All the time, we used the very simple continuous ap- 
proach. Why it is so good? Scaling behavior of networks 
arises from a power law singularity of degree at s = 0, 
i.e., for the oldest nodes, k(s, t) oc s~@. Such, "the oldest 
are the richest", behavior is often perceived as a defect 
of a preferential linking scheme [jjrjf . If one removes such 
singularity or makes it more weak, a growing network 
will be out of the class of scale-free nets. It is the sep- 
aration of power and exponential dependences in P(g, t) 
and p(q, s, t) that makes the continuous approximation 
so efficient. Therefore, we see that scale-free networks 
are quite suitable for the continuous approach. 

Natural boundaries, s — and s — t, are always 
present in growing network. We showed, that even pres- 
ence of strong nodes did not lead to violation of the rule 
"the oldest are the richest" for nodes in the continuous 
part of degree-distribution. Nevertheless, we have found 
a threshold value of the strength of this node above that 
a single node influences evolution of an entire network. 
It captured a finite fraction of all links - giant capture 
- and determines values of exponents, although the net- 
work stays scale-free. These "collective" effects are ex- 
plained by absence of any "interaction distance" in the 
process of network growth. Each node has some chance 
to get a new link. Therefore, in principle, it is possible 
to adjust parameters to direct all new links to a single 
node. 

One should note that main part of our paper has been 
devoted to studying of a one node characteristic of a net- 
work - degree distribution. The same simple character- 
istic is a matter of interest of most of modern experi- 
mental and theoretical studies. This restriction let us to 
apply usually rather general models in which links ap- 
peared between arbitrary nodes: new and old, old and 
old. If we do not study statistics of connections between 
different nodes and global connectivity properties of a 
network, usually the considered problems are equivalent 
to a classical problem of distribution of new particles be- 
tween increasing number of boxes. In fact, in the present 
paper, we studied the question, how a network is self- 
organized into a scale-free structure using some versions 
of stochastic multiplicative processes |2lJ which are most 
known by their "econophysical" applications. Most of 
our results may be described both in terms of the theory 
of evolving networks and econophysics - wealth distribu- 
tion processes. 

Nevertheless, a part of our paper has been devoted to 
problems that cannot be understood using a one node 
characteristic. In this intriguing direction is, perhaps, 
the main aim of statistical physics of growing networks - 
description of topology of evolving networks. 
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new node 



a) 



FIG. 1. Scheme of growth of two standart network mod- 
els considered in the present paper, (a) At each time step, 
a new node with n incoming links is added. Source ends of 
all links are placed anywhere. In addition, m links are dis- 
tributed preferentially among all nodes. This means that a 
target end of each of these links is attached to some node ac- 
cording to a rule of preference, (b) A citation network model. 
n = 0. At each increment of time, a new node with m outgo- 
ing links is added. Target ends of these links are distributed 
preferentially among old nodes. 
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a = 2.0 a = 10.0 

FIG. 2. Change of the structure of the network with aging 
of nodes with increase of the aging exponent a. The aging 
is proportional to r~ a , where r is an age of the site. The 
network grows clockwise starting from the site below on the 
left. Each time one new site with one link is added. 
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Exponential 


Oldest nodes capture 
finite fraction of links 







Scale-free 



A' 



FIG. 5. Phases of the network with the varying strength 
of nodes, G(s) — s~ A . Scale-free networks are realized only 
at A' = 0. 



a 



FIG. 3. Exponent j3 of the mean connectivity vs. aging 
exponent a. Points are obtained from simulations [?]. The 
line is the solution of Eq. The inset shows the analyt- 

ical solution in the range — 5 < a < 1. Note that 1 if 



a — > — co. 
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FIG. 4. Exponent 7 of the connectivity distribution vs. 
aging exponent a. The points show the results of simulations 
[ ?]. The line is the solution of Eq. ( pll| ) with account for Eq. 
(p_2[). The inset depicts the analytical solution in the range 
-5 < a < 1. 
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FIG. 6. Giant capture. Fraction of all links, d/(m + n), 
captured by a single strong node at long times and the scal- 
ing exponents f3 and 7 vs relative strength, g, of the strong 
node. The giant capture occures above the threshold value, 
g c = l//3o — 7o — 1 > 1. Here, f3o and 70 are the corre- 
sponding exponents for this network without the strong node. 
d/(m + n)[g — > 00] — > /3q and f3[g — * 00] — *• 0. 
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Log t Logq 

FIG. 7. Schematic plot of the degree distribution of the 
network with one node which strength exceeds the threshold 
value. The giant peak is due to the strong node. A hump at 
the cut-off of the continuous p art of t he distribution is a trace 
of initial conditions (see Sec. VIII B and Ref. [?]). 
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FIG. 8. Fraction of all links, d/(m + n), captured by the 
component of the network consisting of strong nodes, at long 
times and the scaling exponents /3 and 7 vs relative strength, 
g, of strong nodes. We introduce two sets of exponents for two 
components of the network, ft and 71 - for the component 
consisting of nodes with the unit strength (contains (1 — p)t 
nodes) and ft, and 7 3 - for the component consisting of nodes 
with the strength g (contains pt nodes). Thin lines depict 
the dependences at fixed values of p. Arrows show how these 
curves change when p decreases from 1 to 0. At p — > 0, we ob- 
tain dependences shown in Fig. [| At p — > 1, d/ (m + n) — > 1, 
P g {g < 9c) -> f3 g, P a (g > g c ) -> 1, ft -> Po/g, 71 -> 1+fl/A), 
7s(5 < 3c) 1 + V(A)fl)> 79(5 > 3c) -> 70- 
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FIG. 9. Phase diagram of the network growing in con- 
dition of permanent random damage - permanent deleting 
of random links - in axis 70, — c/m. At each time step, m 
new links are added and — c random links are deleted (see 
the text). 70 is the scaling exponent of the corresponding 
network growing without deleting of links. Curves in the 
plot are lines of constant values of 7. 7 = 00 on the line, 
70 = -2(c/m) - l/(c/m). When 1 + sfl < 70 < 2v% 
there is a maximum of the dependence, 7(— c/m, 7ofixed), at 
-c/m = (1 + ^2)[7o - (1 + \/2)]/a/2. Between the lines, 
—c/m = 7 — 2 and —c/m = (70 — l)/2, the 7 exponent varies 
between the values 2 and 1 correspondingly. The NSD regions 
denotes areas where stationary distributions are absent. 
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